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In conjunction with recent numerical A OqAq + V • vl = "A-gauge" results re- 
ported in a companion paper, we construct an ^ cxd Wilson loop picture of 

T 1 

. A-gaugefixing in which (I)the A-gauge expectation value of a link chain C is the 

O ■■ weighted sum over Wilson loops made by joining to C all selfavoiding chains C dos- 

, ing C. (Il)Weights A~, containing all the A-dependence, are given by the /3 = 

O . ^ ~ 

■ A-gauge expectation value of C. illVjA- equals path-products of coefficients from 

OS ■ 

the trace expansion of the gaugefixing Boltzmann weight. From (II) and (III) we 
\ deduce formulas for /3 = quark matrix elements. We find that Mq"^"^ decreases 

with increasing A; the quark propagator dispersion relation is not covariant when 
' A 7^ 1; and A/ = 1/2 matching coefficients are A-independent. These strong cou- 

pling features are qualitatively consistent with numerical /3 = 5.7 and 6.0 results 
\^ • briefly described here for comparison purposes but mainly presented in a companion 

a 

paper. 
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I. MOTIVATION AND RESULTS 



Traditionally gaugefixing is done only in weak coupling perturbation theory, where it is 
needed to define perturbative quark and gluon propagators. In lattice QCD, gaugefixing 
is unnecessary for computing gauge invariant correlation functions. However, since local 
gauge symmetry cannot break spontaneously the only way to see nonvanishing gauge 
variant correlation functions is by gaugefixing. Accordingly, lattice gaugefixing has drawn 
considerable attention in recent years 0. As described in our companion paper 0, gluon [Q, 
quark |^,|^ and photon propagators effective masses, and wavefunctions have been 
studied in special cases of "A-gauges" 

XdoA° + V ■ A = 0. (1) 

The numerical work has prompted analytical and computational studies of longitudinal PJlO[| 



and topological |TT[ gaugefixing ambiguities and their effects on gluon, quark and photon 
0] correlation functions and operator product expansion coefficients determined from gauge 
covariant matching conditions [p|,p!0|. 

Nonperturbative gaugefixing is a complex subject. Since quarks are confined Mq may 
(or may not) depend on gauge. For example 0, compare the exactly solvable Schwinger 
model in Coulomb gauge to covariant gauges parametrized by gauge parameter ^. (We were 
not able to solve the model in A-gauges.) While the actual situation in dimension D = 4 
QCD and lattice QCD may be arguably different, it is helpful to have a litmus test for 
discarding broad arguments ("Mass is gauge invariant in perturbation theory; hence quark 
mass is gauge invariant.") which do not distinguish between the Schwinger model and other 
gauge theories. In the Schwinger model, quarks are confined but the photon is physical and 
has a mass from the Ua{^) anomaly. We define "effective mass" as the inverse correlation 
length of the zero momentum propagator. While the photon has a gauge invariant effective 
mass (equal to its physical mass), the effective quark mass varies with gauge parameter ^. 
Coulomb gauge — the unitary gauge for the Schwinger model — has no unphysical modes and 
the quark propagator is the physical amplitude for quark propagation. Due to dielectric 
breakdown of the vacuum |]T2[, the quark propagator — the amplitude to have just one quark 



at X 7^ starting with an x = quark — vanishes and the effective Coulomb gauge quark mass 
diverges. In covariant gauges, the presence of unphysical modes (or, alternatively, Gupta- 



Bleuler physical-state conditions on the Hilbert space) ruins this physical interpretation of 
the quark propagator. Covariant gauge quark propagators are not amplitudes for quark 
propagation and do not vanish despite confinement. 

As reported in our companion paper effective quark and gluon masses Mq and Mg 
were evaluated on /? = 5.7 and 6.0 quenched D = 4, color = 3 Wilson lattices in A-gauges 
by matching gaugefixed p = propagators at large t^; to the free particle ansatz 



^lim^Ee^^-^(K.^.^oV^J) = ' (A) J^'''- ' (2) 

where Ej^^^p = 0) = Mj^^^ and Z^'^'^ = 1. The "(A)" superscript anticipates A-dependence, 
although sometimes we will omit it for brevity. The role of background gauge field gauge 
transformations V^; will be explained shortly. The idea of monitoring chiral symmetry with 
quark masses goes back to the Gross- Neveu model [l^, where the flavor Nf — > oo effective 
quark mass Mg = + irrriq + (9(m^) is a continuous, increasing function of bare mass rriq. 
Similarly, matching the vacuum expectation value of the operator product expansion of ipx'^o 
in massless QCD to a free fermion propagator yields for = 3, flavor Nf = 3 and Landau 
gauge fl^ 

Jim M,(,.p^) ~ i%)(l(4)'"(ff„V,)„). (3) 

Numerical fits in lattice QCD to M^^^) = fc^^^M^ + M'^^\ motivated by the CPTh relation 
oc niq between pion and current quark mass f^, yield ~ 2.7(.3) x / MeV and 
M^(^) ~ 350(40)71^61^ at /3 ~ 6.0. At both (3 = 5.7 and 6.0 effective quark and gluon masses 
decrease as A grows [g so that, roughly, /M^'^^ ~ .9 and Mf^'^/M^^'^ ~ .75, all plus or 
minus ~ 15% jackknife errors. 

In this paper, we put forth a A-gauge /3 = 0, color N oo solution of lattice QCD with 
quenched Wilson fermions in an infinite volume lattice. Lattice gaugefixing is implemented 
by a lattice Fadeev-Popov method; we couple the links to a quenched Higgs gaugefixing 
field {Vx} in the fundamental representation of SU{N) 0. Gluons do not propagate at 
/3 = 0, where link fields oscillate randomly. Let us focus momentarily on quark propagators. 
A-gauge quark propagators are Vxipx "meson" propagators in this formulation — as written 
in (0) . An expression for Vxipx meson propagators follows from hopping expanding the Higgs 
and quark fields, doing the (3 = link integrals — which project out all but zero area Wilson 
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loops — and resumming the hopping expansions. However, as described in Section P Q , we 
are unable to resum graphs where Higgs paths recur (Recurrence and other such notions are 
defined in Section [11 A| .) since the hopping expansion weight of such paths differ from same- 
pathlength nonrecurrent graphs. To get around this we resort to a "trace orthogonality 
approximation," which does not differentiate between recurrent and nonrecurrent Higgs 
paths. This approximation is tantamount to taking the N oo limit before resumming 



because, as we show in Section [II Q , there is no difference between recurrent and nonrecurrent 
paths in the brutally truncated — oo limit. Taking N ^ oo before resummation is an 
approximation because the hopping expansion is (apparently) not absolutely convergent at 
infinite A^. 

In the trace orthogonality approximation, only selfavoiding quark paths contribute to the 
P = quark propagator because the {Vx}, being quenched, dress only quark paths without 
internal loops. Hence nonselfavoiding quark paths are suppressed by infinite string tension 
at /5 = 0. However, since are unable to (re) sum over only selfavoiding quark paths, we make 
an additional approximation and sum over all (for technical reasons) nonbacktracking quark 
paths. When we do this we find, as shown in Section |T|, that the (3 = 0,r = 0,N^oo zero 
momentum quark propagator pole Mjj'^^ is analytic and linear in m„ as m„ 0. Expanding 



yields 



/Tyr^ — T ril + 14D-25A2 A<i; 

(,)^ V2D-1 -2' 

9r.^9n-7U I -25+200A+224(D-l)A2 ^ ^ 1. ^ ' 

^ ' \ 4{4A-1) ^ - 2' 



499-1474D+1600£>^+(1376D-751)A^ 
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I -751+1376D+6008A-11008DA+(-4032-1568g+25600i:'^)A^ X > 1 
[ 8{32D-7)2A(4A-1) ^ — 2" 

Note that the A-dependence of M^f^^ is qualitatively the same whether D ^ oo (where 
the number of recurrences and selfintersections are negligible) or D = 4. This suggests 
that recurrences and selfintersections are not responsible for A-dependence. Therefore we 
generally quote formulas for all D when we might be safer with their D ^ (x limits. In any 
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case, we do not find any qualitative difference between tfie A-dependence of M^^^ at finite 
and infinite D in our approximation scfieme. 
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FIGURES 

FIG. 1. A-dependence comparison of numerical quark masses to the /? = 0, — > oo For- 



mula (4b) evaluated at D = 4. The numerical masses are rescaled so that at A = 1 all data points 



are normalized to Mc^\ Ref. Q provides details of the numerical simulation. 

As depicted in Fig. |I], Eq. (^ ) mimics the A-dependence of the (3 = 5.7 and 6.0 numerical 
data. Most of the change occurs between A = 1 and A = 2, and Mj^^^ stabilizes to a 
nonzero value as A ^ oo. This qualitative agreement between strong coupling and numerical 
behavior helps give confidence that the numerical A-dependence of Mg is not a finite volume 
artifact or due to details of how Mg is extracted in the numerical simulations. 

For technical reasons we compute the quark propagator dispersion relation not at r = 
but at r = 1, where it is 

= Mj^'^' + Y.t f^ + 0{a') , ) ^ ±tpo , (5a) 

i=l 

, (5b) 

Since g^^X 7^ 1) 7^ 1, the quark dispersion relation is not covariant and, hence, the propagator 
is not free particle-like except in Landau gauge, drops out if p = 0, where the quark 
propagator is indistinguishable from the free particle propagator. (However M^^^ remains 
A-dependent.) 



As described in Section [HI Q , matching coefficients for sd subtraction of A/ = 1/2 Rule 
operators in the f3 = 0, N ^ 00 limit assuming trace orthogonality are given by 

«o±[rir2] = ((^riFa.ss T ^g/l^y /rfra) (^0^0) (6) 



where Fierz coefficients /rfr2' given in Eq. ( 87b ). Since is gauge invariant these 



/? = matching coefficients, which are directly related to physical continuum decay rates, 
are gauge invariant. At /? = 6.0, ao±[rir2] also seem to be A-independent up to statistical 
errors W\- 
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II. LATTICE GAUGEFIXING 



As exemplified by several models — notably QEDa+i — the thermodynamical limit of 
strong coupling lattice gauge theories may correspond to different field theories than the 
(3 —>■ oo ones. In particular, proving gauge dependence of quark mass in the lattice Schwinger 
model would not reveal much about quark mass in the continuum Schwinger model since 
the lattice Schwinger model has a phase transition in (3 |T^. (This transition, if the critical 
point is unique, doesn't ruin confinement since Schwinger model quarks are confined at both 
weak and strong coupling.) Analogously, the following P = 0, N ^ oo solution should be 
viewed as a toy model and not something necessarily related to QCD. 

The A-gauge lattice expectation value of a lattice operator 9 is 



(6) 



(7a) 



where 



[de']e 



6 e {f,v,u}, 



(7b) 



x,y 



(7c) 



S-[n] = -Y,^-^ trfi, „ , n e ±{0, ■■■,D-1}, 



J-n = Jn = ^n/ii ^fj. = ^i ' ' ' ^ ^) ) ^x,n = VxUx^nVx+h^ 



(7d) 
(7e) 



5"[f/] = = /3 ReJ^trD 



(7f) 



The sum in (|7|) ranges over all lattice plaquettes trD. In ( [7a| ) quarks are quenched by 
choice. Consistency with the Fadeev-Popov method requires {Vx} to be quenched — inverse 
partition function z^^ = z:^^[U] plays the role of lattice Fadeev-Popov determinant. Since as 



described in Section |TT| the hopping expansion turns [0]/ into a composite operator of links 
and transformations, in this Section we focus on the gaugefixed expectation value (VyyAwV^) 
of continuous link chains yA^. The strategy is to integrate out {Vx} to get an expression for 



(8) 



This integration is simple enough to be practical because gauge symmetries relate V to 
continuous link chains and suppresses all disconnected link chain loops. Then gaugefixing 
can be viewed as an operator insertion of V into 



{Vyiy^wW^) aj - [Vyy^wVlU - iy^w) l3i ^ al3;ij[U] 



(9) 



the usual gauge invariant lattice expectation value. 
We will show that in the ^ oc and N ^ oo limit 



(10) 



where the sum ranges over a complete set of selfavoiding continuous link chains Or- 
thonormality of the ^Ay with respect to inner product 



(GilGa) = J [dU] ^ tr(ele2) , [dU] = I[dU,,, 
identifies coefficients A ~ with 6 = A-gauge expectation values: 



(11) 



(12) 



If y = w = X, the only continuous selfavoiding chain is = 1- Since V^V^ = 1, Eq. ([I2| 
implies = 1/N. Hence [(yx)ai3{Vj)ij]y is independent of [U] and 

Eq. (|13D implies gaugefixing does not affect closed link loops Q)^ since 

1 



(13) 



(14) 



A. Link Operator Definitions 

A continuous link chain or randomwalk path may selfintersect{tonch itself at right an- 
gles), recur (touch itself at 0° or 180°), or feacA;tracA; (immediate recurrence as in the sequence 
"■ ■ ■ Ux,^Ux+ji-^i ■■■")• A selfavoiding link chain or randomwalk path does not selfintersect 
or recur. 
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A tree T is a continuous randomwalk path which is nonbacktracking but may be self- 
intersecting and otherwise recurrent. A branch i? is a continuous, possibly selfintersecting 
and definitely recurring randomwalk loop enclosing a zero area minimal surface. A tip of 
i? is a backtracking subsegment of B. Every B has at least one tip, perhaps more. Every 
randomwalk path is a sequence of trees and branches. 

Continuous link chains extending from y to w, not necessarily straight or selfavoiding, 
are denoted by 

yA-u) = Uy^^UyJ^fi^u " " " U^_g^^ , jy^s]^ = ^Ay (l5) 

where "• • •" denotes a continuous but not necessarily selfavoiding link chain. Selfavoiding 
link chains yA^ are accented with Examples of recurrent and selfintersecting link chains 
are □., a continuous chain which traces out a unit box and recurs at one side, and ^, which 
wraps around a unit box and selfintersects at one corner. Link chain loops ^Ax are also 
denoted by Q- The unit square chain is Q = whose trace is the plaquette. If Q begins 
and ends at x, then it is denoted Qa;- 

If the path traced out by Q^- is a branch, then it is denoted Of - Since Ux,^Ux+ji-^ = 1 

Of = 1 vof • (16) 

Factorizing link chains into a product of tree and branch subsegments yields 

yAw — Oj/ y^x Ox ■ ■ ■ Oz ^^luOto ~ y^w ■ (l''') 

Hence we will assume that yAw = yA^ unless specifically noted. 
By SU{N) identities such as Schur's lemma, 

J dU D^iU) D%\U) = ^r^f'"' S^pS,^ , (18) 

continuous link chains are orthonormal with respect to (|11]). However some disconnected 
chains mix with each other. For example, a nonorthonormal basis for gauge invariant link 
operators is 

S, = {l,{tr0},{tr0tr0'},---}1, trl = AT (19) 

where {0} is the set of all nonbacktracking link loops in the lattice. It is possible to 
orthonormalize Bo by constructing irreducible combinations out of its link operators which 
mix. 
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B. Residual Gauge Symmetries and the Link Expansion 

QCD gauge transformations are 

While S"" breaks QCD gauge symmetry, S*", S-^ , and are invariant under 

U,,^ ^ R^U^^^rI^^, V, V,Rl tP, ^ R^ij, . (21) 

We will refer to (^Dj) as "Q" transformations and ( PT| ) as "i?" transformations, i? is a 
symmetry of [Q]e\e<^{f,v,u} (and hence (6)) whereas Q is a symmetry of only [Q]e\ee{f,u}- 
Elitzur's theorem [ij any /2-variant operator such as and [4^^ has zero (0) expectation 
value, /^-invariant operators, including Q- variant ones like Qx,tJ., are not suppressed. 
Global transformation "L^y" where 

Ux,f, , ^ ^xyK , tpx^tpx (22) 

is a symmetry of ['S>]0\ee{f,v,u}- Lxy is equivalent to global color symmetry. The XY des- 
ignation is because in the (3 ^ oo limit the system approaches the XY model, which has 
Lxy symmetry. 

Under Ux,^ — >■ U^^^ = QxUx^fj.Q^^x+fi partition function z.u\U] obeys 

z^[U'] = J [dyje-^"!^'] = / M(Vgt)]e-^^Pl = .^p]. (23) 

As a Q invariant functional of links Zy, the inverse of the Fadeev-Popov determinant, is 
expandable in terms of gauge invariant link structures. Such structures are comprised of 
links joined together into closed networks with gauge invariant bonds 6ij and ei-^...ij^. Let 
us take N oo or U{1) and throw out the latter. Then following (0), with O 7^ 1 ^^'^ 

zM = i^ + T.Z?tTO+ Y: ZpO't,QtrO' + ...)z:, (24a) 
^ / [dU]zy[U] = ZI+Y: ZpO\, o tr + ■ ■ • (24b) 
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The reason for factorizing out Z° in (|24a|) and the reason it can be expanded in terms of 



zero area (colorless) link structures is because Z° is the colorless part of Zi,[f/]. Since it is 
colorless, Z° survives at /? = 0. 

If W is in the fundamental representation and 'y^^'^ the character of irreducible represen- 
tation z/, the character expansion of e^^*^*'''^ is given by 

h^{z) = K{z) = [dW x^'HW) eiMH/+H/t)^ ^25a) 



e-Mw-^^^) = J2bA,)^M^W), fu^^. (25b) 



Since by (|7ie-^^'[^l = n.,^e 



^-s^M = zi Y: IlNU.JJ,)x'^''-^\^.,,) (25c) 
where the reason for factoring out Z° will be apparent shortly and where 

fdW O pi2{tr(W+H/t) 

1|/„W1|<1, lim/„(;) = -i-, Mc^) = l (finite JV). (25e) 



The first inequality in ( p5e| ) follows from ||triy || < N; the second from expanding e^^''*'"'^ in z; 
the third since the integrands are dominated by MaxRetrW^ when z ^ oo. Integrating both 
sides of ( p4a| ) over [dU], interchanging integration order — > [(iy][(if/] and changing 
variables [dU] to [dQ] yields 

= I [dU][dV] Jle^-^^*--^-- = X{h^{J,) . (26) 

x,n x,iJ, 

LxY symmetry implies 

[{Vy)afsiV^hl = 6^,U^[U;w,y] . (27) 

Since t[U;w,y] transforms under local R like [(K/)a/3(Kl)jj]t)5 at y it transforms like N, at 
w like N, and the links away from y and w must be bonded gauge invariantly either by the 
SU (N) identity tensor or the completely antisymmetric e tensor. 

To appreciate that (^) is not a selfevident result — it's valid only in the N ^ oo limit 
and orthogonal trace approximation to be described — consider the U{1) case where 
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cos 



X "x+fi I H^x,fj,J 



X,fJ, 



The U{1) character expansion is given by 

oo 



-ind 



= {{nx,t,}\^n^,m = 0} , 



x,m ^x,y 



^x,w } 



The U{1) identity used to obtain (|28d| - p8g| ) is 



m, rr^i m. 



where x^'-^'H^x.^) = X^'^^-^-'^H^La.m)- 

Because character coefficients do not satisfy 

u{z) = {u{z)T , 



(28a) 

(28b) 

(28c) 
(28d) 

(28e) 
(28f) 

(28g) 
(29) 



(30) 



hnk loops in the numerator of [K/^^Ji; which touch open chains ^Ay are not cancelled by 
corresponding z^ loops and (0) does not apply to U{1). Formula (|TUp only applies to 
systems consistent with both (^) and (|30|) . 

The SU (N) version of ( ^9l) is complicated because there is more than one way of making 
a singlet out of i'-d+i CS) ■ ■ • C?) J^d-i- SU (oo) traces, on the other hand, obey an analog of 



my- 



D-l 



fdV. n tr"--(\40.,„)tr'--(0t^„\/i) 

m=-D+l 

= ^(^x,D-l — L,D-l) 



(31a) 



^x^S- ^ ^ jx=0 



^ ^^g-[tr(0.,M.o-)Oi,^^(, 



N 
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where 



^ X — UM ' X — U 



m 
m 



/i> 0; 
-V < 0; 



(31b) 



m = —u: 




(31c) 



(l-D < J < n^, 



n'X,m — ^ ^ ^x,m'j 
m'=-D+l 



l-D', 



(31d) 



D - 1 r2^,D-2 < j < n^^D-i- 

Sp is the permutation group on p elements. Eq. ( pia| ) says that the numerators and denom- 
inators of t[U ; w, y] are comprised of chains made by joining hnks from the trace expansion 
of e""^" in all possible permutations consistent with {rix^m, ix,m} conservation at each site. 
In Section |11 C| we will find that relevant coefficients of the link expansion satisfy ( |30D and 
that the SU (oc) analog of ( p8cl[ - p8g| ) is consistent with Eq. (|T0|). 



C. Trace Expansion and Orthogonality Approximation 



We will call the use of (|31a]) the "trace orthogonality approximation." The reason its 
use is an approximation even in the N —>■ oo limit is as follows. If 2 < < cxd, an exact 
group integral formula is 

1 



"2/^2 _ I 

{^hl3i^jiai^i2f32^j20i2 ~ 'T7^iil3i^jia2^i2l32^j2ai) "I" ( 



«! a2 



(32) 



Since the second and fourth terms in the RHS of ( ^2]) are suppressed by 1/A^ relative to the 
other terms, (p^ ) agrees with ( pia| ) in the N oo limit. However, using (|3^ ) in place of 
( piap to compute the leading 0{J^) hopping expansion term of Au, the Eq. (|iy) coefficient 
where is a continuous five-link chain which traces out a unit square and recurs at one 
side, yields 



1-1-1 



(33) 
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The first two terms of tfie middle expression come from the numerator of t[f/; w, y]; the third 
from a denominator plaquette and a numerator link. Since J oc N (as described below), 
Aa = 0{1/N) if we take N ^ oo after resummation of the trace expansion. On the other 
hand, An = if one uses Eq. ( ^la| ), vahd to leading 0{1/N), to integrate out the [dV] (as 
described below). The latter is tantamount to taking ^ oo before resummation. The 
discrepancy arises because the leading 0{1/N) contribution to Aa comes from subleading 
terms of ( pla| ) which have been thrown out. We have not been able to improve the trace 
orthogonality approximation to account properly for coefficients of recurrent chains. 
The orthogonal trace approximation leads to the following for 

Z[L, K]= JdW eh^r[LW+W^K]^ z^j^ J] ^ ^34g^) 

Expanding in L and taking N ^ oo, imposing 

J dW tr"(Liy)tr'(iytir) = n\ [^^^^^]" Sni {N oo) (34b) 
and resumming gives 



;=o n=0 



( p4b|) is a special case of (|31a|) . We change notation from Z ^ Z to emphasize that the 



interchange of integration and resummation leads to a discrepancy, described below, between 
Z and the original Z. 



Using Z as a generating function yields ("/ /" is defined in ( |25d|) . 



F^d.:H,r = /Win ■ ■ ■ Wl^Wk^o, ■ ■ ■ Wt^oJ , (35a) 



2^ =''-w^ =e'=fr=%/«(^). (35b) 

OL^i ,L=j oKji ,L=j 

\k=j ik=j 



fNi-J) = t4 « 2iV' N oo). (35c) 



2N 

Since following ( p4c|) 



(9"+' log Z 



dK^^i, ■ ■ ■ dKj^i^dLo^ki ■■■dL 



(if n or / > 1), (35d) 
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(35e) 



by induction. Combining ( |35b|) and ( |35e|) yields 
1 



/tr"iy tr^Wy = fp-^J) (n > 0, / > 0) 



(35f) 



Eq. ( P51|) is our preliminary SU{oo) analog of (30) which we will now improve. 

Formula ( |35cj ) violates the inequality of ( |25e| ) if J > 2N because Z is not valid to leading 
0{N). Rather, following ( |34b| ), it is only asymptotically valid as ^ oo. Since (|34c|) is 
analytic in tr(Li^'), the violation implies Z ^ Z. This discrepancy arises from applying 
to (|34c|) on n — s> oo (including n > N) terms before resummation. 



An exactly solvable limit is ^ oo and J oc N. In this double limit, achieved by 
identifying the gauge parameter ^ with ^ oc 1/A^ so that 



AT ^ oo, C 



^ 

' 2N 



X 



(36) 



the A^ dependence cancels out leaving a nontrivial A-dependent solution. The exact gener- 



ating function is [|T^] 



(37a) 



Ik 



-'^ — 2' 



- \ log( 



'iK > 



2' 



(37b) 



lim /jv(J) 



5^0 



dw 

d£ i=K=\ 



(37c) 



By induction 



i,/tr^l^/-/.(.)^i,|^^0(i,). 



(37d) 



(37e) 



Since the LHS terms in ( |37d| ) and ( |37e| ) are 0(1), property ( |351D is valid when J = 2A^A in 
the N ^ oo limit. 
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The link expansion of t\U\w, y] follows from trace expansion 

tr"fi tr'^]"^ 



oo oo 



e^^"^^ = EE-K^;-)^^- (38) 

n=0 «=0 '■• 

Taylor expansion of e^^'^*''^, which defines the exponential operator, implies r(n, /; J) = 
(J/2)""'"^ In fact, this identification is only asymptotically consistent with using ( pia| ) in 
the limit J « 2N and A^^oo. AsJocA^^oo, integral formula (|31aD cannot be 



interchanged with resummation. If we insist on applying (31a) before resumming the trace 



expansion, r(n, /; J) must be redefined so that correlation functions evaluated using (|31a]) 



and (|38D are consistent with correlation functions evaluated directly from Z of Eq. (|37ai) 



The matching condition 

t "'"'"l'--'^ = h{J) iV fUJ) = t "'"-'''"•^^ (39) 

n=0 n=0 

implies when N ^ oo 

(I J « 2N; 

r(n,Z;J) = iV'+"/^+"(J) , , (40) 

The 5f/(oo) analog of (^, 

r(n, / + A;; J) = r(n + k, l; J) = f%{J) r(n, /; J) , (41) 

will be referred to as "factorization." 

Applying Eqs. (|31a|) and ( pSD leads to the SU {oo) version of (P8d| - |28g]) . The [dU] integral 



of Eq. (O) for A ~ projects the the numerator and denominator of t\U] w,y\ down to the 
colorless part except along w^l^. At other sites the denominator contributions cancel the 
corresponding numerator contributions. Each link of coming from a (|31a|) integration 



costs a factor of 1/A^. Additionally there is an overall factor of iT{w^\w^y) / N"^ = 'i-fN. 
Hence 

^ i = ^TT^ ru I J) =Tr n fNi-Jm) 42a 



where 



( 1 if {x,m} e uj^y] 

A{n,x,m} =1111' P^^rn = { ^ (42b) 

{nx,m.} {x,m} [0 if {x, 777,} ^ w^y 
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Let a(0) be the coefficient of link chain in the numerator or denominator of t[U; w, y]. 
Consider again An where is a hnk chain which traces out a unit square and recurs at one 
side. The chain □. does not contribute to An in the same manner that ^Ay contributes to 
A jjj because, due to factorization, a(n.) = (^{~) ^■(n) so that 

q(-) + q(D) l + a{n) 

; ; — ^ — = a — ) ; — - = a — ) . (43) 

l + a(n) ^ ^ ! + «(□) ^ ^ 

(^) says that the □. contribution to An has aheady been counted in A^. This imphes that 

recurrent or ^= . ^44^ 

sclfintcrsccting chain 



As described in the beginning of this subsection p^), the trace orthogonahty approximation 
result, is inconsistent with the leading expansion contribution. The problem is that the 
leading contribution to An is not from the leading term of (|3^). In ( |33D the middle "—1" 
comes from the second term of (^) — a term which is neglected in (|31a)) . Hence it does not 
appear in An if ( piaj ) is applied to do the [dV] integrals. 



D. Wilson Loop Picture of Gaugefixing 



Eqs. (^), (|T0|), and ( 42a ) imply that as — *• cxd within the trace orthogonality approxi- 
mation 

= ^^AVyy*^.V%=, = n ^ fN{J,) , (46) 



where /at is given in ( |35c| ) if is finite and in ( p7c| ) if ^ 1/N. These formulas have the 
following interpretation: 



By Eq. (]45| ) the gaugefixed expectation value of a link chain yA^ is the weighted sum 
over all Wilson loops made by joining to a self avoiding link segment ^Ay. This 
holds at all /3. 

By (5B|), weights A ~ are proportional to the 3 = gaugefixed expectation value of 
wAy. All the gauge dependence is in these weights. 
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• A ~ is a path-product of e trace expansion coefficients. In the orthogonal trace 
approximation ( |34b| ), nonselfavoiding operators do not contribute to and the 
/? = expectation value of nonselfavoiding links vanish. 

• Following Eq. ( p!] ) P = A-gauge expectation values satisfy factorization. If Psa is a 
selfavoiding path 

( n ^y,n)\f5=o= n {^y,n)\p=o- (47) 

{y,n}ePsA {y,n}^PsA 

These results permit us to read off the Wilson line A-gauge expectation value, which is 
proportional to the heavy quark propagator §]. Asymptotically it can be parametrized as 

J™. ^t^(^(*-^^^^,.^+t.8^(L^) = ^^e-^'^'*-', X ^ {tE, 0). (48) 
V^^'* can be interpreted as V^^^ = M^^g^^^y — rrig. At (3 = and N ^ oo, 

Zh = 1, 4'^ = -log[/^(Jo)]. (49) 
Hence V^^^ decreases with increasing A. 



III. HOPPING EXPANSION AND RESUMMATION 

In the N oo hopping expansion of [0]/, the ipx operator hops from site to site leaving 
behind a trail of 

(r - fL) KBUy,n , riG {±0,---,±(D-l)} (50) 

factors to mark its randomwalk path P. Pauli Exclusion is irrelevant because quark quench- 
ing effectively requires quark paths to selfintersect indefinitely.S 



^The author thanks C. Bernard for pointing this out. 
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Therefore the quenched background field quark propagator is 

[^.V^o]/ = E n {r-ii)KBUy,n (51) 

P {y,n}eP 

where {y, n} are the locations and directions along P. (Assume implicit pathordering when 
appropriate.) Following Eq. (|13), decompose the sum over P into a sum over trees and 
branches. The gaugefixed expectation value of branches is trivial by (0). /3 = expectation 
values of nonselfavoiding trees, as discussed in Section |II Q , are suppressed in the orthogonal 
trace approximation. Thus following (|4^) Eq. (^) is equivalent to 

n = (V;^.?^o^J) = E [ n fN{Jn)] E W {r-mB (52) 
{Tsa} [vMe-TsA {By\y&TsA} {yM<^P 

where P is the selfavoiding tree Tsa going through sites y with a branch By at each site. 

As there is no technique for summing selfavoiding trees as specified in (^) , we shall sum 

all trees T. This approximation is justified if D > 4 by the fact that the quark mass pole 

is determined in the t > L ^ oo limit. In this limit the average number of selfintersections 



and recurrences 17 



s = L X (path density) ~ j^— = L^"^/^ (53) 

vanishes if D > 4 for randomwalk paths. (Nonbacktracking randomwalks are presumably 
at least as selfavoiding.) This suggests that the number of nonselfavoiding paths, being 
negligible, may not affect the quark mass pole. 

Whether selfavoidance affects the propagator pole or not is analogous to whether 

oo 

/(a, x) = aL , lim ul = (54) 

L=0 

guarantees f{a,x) is of the form f{a,x) oc + finite. A counterexample satisfying (|54D 
is = a"^ + 1/L in which case a second pole exists at x = 1. Hence anomalous quark mass 
poles or mass shifts stemming from the inclusion of nonselfavoiding paths is a possibility 
which cannot be ruled out. 

There are three effects we wish to investigate: (i) spontaneous chiral symmetry breaking, 
(ii)gauge dependence of the dispersion relation, and (iii) proximity of the effective quark 
propagator to free particle form. Since the source of (i) is certainly different from (ii) 
and (iii), we may pursue them separately. If the dispersion relation is gauge dependent 
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at {r,f3), there is no reason — barring a critical point — to believe this property would not 
persist to other (r, /5) values since (r, /3) are unrelated to gaugefixing. To study (ii) and (iii) 
we set (r = 1, /3 = 0) where we can solve for the full quark propagator and its dispersion 
relation. In this limit, spontaneous symmetry breaking is absent. To study the combination 
of spontaneous symmetry breaking and gauge dependence, we set (r = 0, /? = 0) where we 
can solve for M^^^^ 



A. — > oo, r = I, P = Quark Propagator 

Since (1 — 7^)(1 + 7;^) = 0, r = 1 quarks cannot backtrack and the hopping expansion 
cannot generate any zero area Wilson loops. Hence at /? = 

(?^o^o)^°"* = -2CKb , C ^ 2^/2iV , Ks ^ ^, ^ ^ , (55) 

2[mq + D) 



where ■?/)^°'^* = y/2KBip and C is from the Dirac and color traces. The r = 1 pion mass 
obeys ll 

,,,,, , l-AKl(D + l) + 16DKi 

= 1 + 8A1(1 - 2(D- l)Kl) ' ^""^ 

which has = l/A at D = A. 

Upon the approximation Tsa —>■ T in Eq. (|5^), the sum over trees is implemented by 
solving the recursion relation 

P,{x) = -Kb 5.,o - KRY.fN{Jn) mx - n) . (57) 



Fourier transform Pq{x) oc J2pG^^^Pq{p) produces 



PM = -- 

90 



[mR + 9^,(1- cosp^)] -iE^lf^Qf, sinp^ 



(58) 



For z G {1, ■ ■ ■ , il* — 1} is given in Eq. (|5b| ). M^^\ the p = pole of Pq{p), is related to 
continuum pole rriR by 

sinh M^^) = ruR . (60) 
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In Landau gauge where = 1, Pg reduces to free particle form with uir = niq/go. 
Otherwise, Pq obeys a noncovariant Dirac equation which in continuum language is 

(*E^.7A + ^fV.(a:)=0. (61) 

Eq. ( |6lD corresponds to dispersion relation (|5aD . 

When K — >■ = 1/4, M^'^) is negative — for nonnegative values of M^'^^ there is no 
massless pion. This reflects the absence of chiral symmetry and spontaneous chiral symmetry 
breaking at (r = l,/5 = 0). Thus we turn to (r = 0,/? = 0) in order to examine the 
A-dependence of Mj^'^^ in the presence of spontaneous chiral symmetry breaking. 



B. iV^oo, r = /3 = Quark Mass 



When r = 0, by virtue of 7„ 7_„ = 1, the Dirac matrix of each branch is 1 and the sum 



over branches in (^2]) is equivalent to renormalizing Kb to [|19 



Kn = KbW{Kb) 



(62a) 



W{Kb)-, the weighted sum over all branch configurations at one site ||19| , |20| , is deduced as 
follows. The number I{L) of length L > 1 "irreducible" branches, branches with a single 
base stem, is 

L-l 



1{L) = (2D-\)Y. Y: m^TT^Hh) (/(0) = 0), 



2D 



the sum over all arrangements of its irreducible subbranches. At any site 



(62b) 



w{Kb) = i + J:{-kI)'^i{l) = — 



L=l 



where 



Wi[Xj 



r — 1 -L n ; 



L=l 



2^W/(X) 



(62c) 



(62d) 



Definition (|623|) factorizes the RHS of (|2§) to give the RHS of {^§2^. Equating the RHS 
of (|62d|) to the LHS gives a polynomial whose solution yields 



l/{mq + ^ml + {2D-l)). 



(62e) 
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Eq. (|2|) reduces (|5|) to 



n(^) = E n {-KR)ilfN{Jn) 
TsA {njeTsA 



(63) 



In this approach gauge invariant (r = 0, /3 = 0) meson propagators are |1T9| , |21 



T {n}eT {mjeTt 



(64) 



where T may be nonselfavoiding since no gaugefixing is involved. Since trees are 
nonbacktracking — a backtrack makes a branch — its Fourier transform Dm{p) obeys non- 
backtracking randomwalk recursion relation 

ii[D^T] = C-KlY, tr[7„r7„(l - 5„,_^)D„e-*P™] (65) 

n 

where F = 75 for the pion and 73 for the rho. The commutation of F with the 7^ splits 



from 



Ml 



niq + 0{rn? 



(66) 



(67) 



While Mp does not vanish as ruq 0, M^^ vanishes and the pion (actually 2 pions) plays 
the role of Goldstone boson for chiral symmetry breaking. The same sort of arguments 
give |2gj21| 



(tAoV^o) = -C- 



DJ2D - 1 + m2 - (D - l)ml 



C = 2^/2;v 



(68) 



for the chiral order parameter, which is nonzero when niq — >■ 0. 



1. Nonbacktracking Approximation 

Since including backtracking overcounts branches, already summed by Kr, we will en- 
force nonbacktracking but otherwise permit recurrence and selfintersection. The sum over 
nonbacktracking trees is implemented by 

P,{x)n = -^^^ - Kr fNiJn) f^T.i'^ - 5™,-n)n(x - u)^ (69) 
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where Pq{x) = J2nPqi^)n- Subscript n indicates the direction from which site x was ap- 
proached. The Fourier transform of Eq. (KOf) produces 



PMn = + T.MnmiE,p)P,ipU , (70) 



m 



Mnm[E,p\ = -Kn fr,{Jn) fl{l- 6m,-n) 6*^" , (71) 

where is a 2D x 2'^/^] matrix. 
Effective quark mass Mj^^^ obeys 

det(l - A^[Mf\0]) = 0. (72) 

Exphcit solution by Mathematica© aX D = 2 and D = A reveals that, despite being a 2D 
polynomial in e~*^« Eq. (|7^) has only two distinct solutions corresponding to iM^^^^ The 
positive energy solution reduces to 

" ' ~ 2(1 +g^K-)goKn ^^^^ 

where g = gi foi i E {1, D — 1} . In the absence of gaugefixing g^ ^ and M^'^^ oo. 
In the absence of renormalization and perfect gaugefixing, — > l/(2mq), g^ 1 and the 
free particle relation sinh Mj^^^ = nig is recovered. 

Following Eq. (^a]), the chiral limit is extracted by replacing Kn with rrig according to 
Eq. ( |62e| ) and expanding about rriq = 0. li d = 2D — 1, 

^(A) ^ 1 + ADjD - 1) - (3 + ADjD - 2))g' - d{gl + g^gl) 

Vd{2d + 2r)go) 

^(A) ^ 1 + ADjD - 1) + di2Dg' + g^) - jd - 2)g' + jg' + 3d - l)g'gl 

2{d + ry9o ■ ^ ^ 

The choice Aj = A = 1 (for i = 1, - ■ ■ ,D - 1) yields (|^) and (|c|). While M^(^) is continuous, 

its second derivative is discontinuous at A = 1/2. 

2. D ^ oo Limit 

In the D ^ oo limit, there is no difference between selfavoiding and random paths. In 
this limit, the leading 0{1/D), r = effective quark mass is straightforwardly expressible 
in terms of the chiral symmetry parameter: 
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Mf) = ^ , = (75) 

l6(?/'o^/'o) go 325(0 

l/gQ^X) = l//Ar(2A^A) is a continuous, monotonically decreasing function whose 2"^^ deriva- 
tive is discontinuous at A = 1/2. 



C. Matching Coefficients 

CaU the weighted sum of (naive) lattice operators whose matrix elements reproduce 
matrix elements of a continuum QCD operator the "lattice representation" of said continuum 
operator and the weights "lattice matching coefficients." Verifying the AJ = 1/2 Rule in 
lattice gauge theory is a longstanding unsolved problem because it has not been possible 
to determine all the lattice matching coefficients of the continuum electroweak Hamiltonian 
responsible for K ^ n matrix elements, related to — > tt+tt^ matrix elements by chiral 
perturbation theory. Specifically, we are interested in operators |^ 

0±™°*[LL] = 4/^" + J2 zl'^'0i}^''[TiT2] , (76a) 



li 2 



where 



o^nrir^] = Fr^V. ± i^i^^., (76b) 
-^riFa = ^orifio ■ gor2go, -f^FiFa = ^origo ■ qo^2do (76c) 
T1T2 G {LL, SS, PP, TT, LR} , J(0) = s(0)rf(0) . (76d) 



The RHS of ( [76a| ) is the combination of lattice operators required to reproduce K —>■ ir matrix 
elements of continuum operator 0±™'^*[LL]. "Matching" coefficients 2;^^^^ are determined by 
gauge invariant numerical methods. The problem at hand is determining z^. Proportional to 
(ms + mrf)a~^, the z"^ coefficients are beyond weak coupling perturbation theory or the usual 
lattice methods. As described in Ref. f^, we use lattice gaugefixing as a technical device to 
determine values of these matching coefficients. The idea is to nonperturbatively replicate 
what is done in WCPTh, that is, to make the lattice representations of the continuum 
operators reproduce an ansatz, motivated by flavor symmetry considerations, for continuum 
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quark correlation functions. This gauge covariant matching condition imposes constraints 
on the lattice matching coefficients sufficient to determine them. We find that matching 
coefficients numerically determined in this way are A-independent (within jackknife errors) — 
as they must be if they ultimately contribute to the matrix elements of gauge invariant 
continuum operators. 

The matching condition invoked in Ref. ^ 

(s|0±™"*[LL]|rf) = , (77) 

the parity even quark-equivalent of the parity odd hadronic renormalization condition that 
K vac, implies that 

4 = - H ao±[rir2F£'^' > (78a) 

{s\0±'^''[T,T2]\d) = ao^[r,r,]{s\f^'\0)\d). (78b) 
Since quarks are not part of the physical S'-matrix, Eq. (|78b| ) may give a value of ao±[rir2] 



different from its physical sector value. In the Schwinger model, |]T0[ not all matching condi- 
tions can be satisfied. For the satisfiable ones, different gauge variant matching conditions 
(relations between different gauge variant correlation functions) may lead to different (or 
similar) values for matching coefficients. The differences stem from unphysical gauge vari- 
ant modes due to gaugefixing ambiguities — gauge variant operators transform differently 
under residual gauge symmetries. Nonetheless matching coefficients are ^ independent and, 
thus, plausibly gauge invariant. Matching with different gauge variant correlation func- 
tions correspond to adopting physically inequivalent definitions of the matching coefficients. 
Hence the use of quark matrix elements (as opposed to, for example, diquark matrix ele- 
ments) must be justified in a physical way before proceeding — as done above Eq. (|78a| ) for 
this AI = 1/2 Rule example. 

We have two purposes in this Section: (a)to understand how A-independence of tto±[rir2] 
comes about; and (b)to derive a consequence of /3 = factorization which can be compared 
to /3 = 5.7 numerical results. This comparison lends insight into why (3 = formulas mimic 
so closely the /? = 5.7 and 6.0 data. 

Following ( [78 bl ) we are interested in A-gauge quark correlation functions of 0-|-^'^"[rir2] 
and J''^**(0). Because of (3 = factorization, quark propagators in the fermionic Wick 
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expansion of correlation functions do not interfere with each other. Hence the /3 = corre- 
lation functions are given by replacing the background field quark propagators in their Wick 
expansion with the /3 = value of (Vxipxi^oVQ). 

Let {a,b} be color and be Dirac indices and define 



X 

where Q G {-D, S", C} corresponding to ■?/' G {(i, s, c}. If ty = —t, by (|6 



(79) 



(80) 



By R symmetry and 



(81) 



In fact, a stronger statement can be made for x- Let "/3 = graph" refer to any single 
hopping expansion graph of [V'oV'o]/ which contributes nontrivially to x P = 0. Then, 
since whenever encloses zero area Y[{y,n}€0^ ^y^^ fi = I and the [dU] integral in ( ^1]) is 
trivial, 



[(^0)^(^0)^]/|,.0 graph ^ ^^'^'^^ ■ 



(82) 



This color diagonality of [ipo'^'olf permits factorization of correlation functions containing 
[i^Qipolf "bubble" contractions at /3 = 0. While [ipoi^olf graphs are not diagonal at /3 > 0, at 
f3 = 5.7 



[(V'o)a(^o)fl]/ « ^"^^^(3 + 0.1% fluctuations. 



On a typical lQby24 gauge configuration with X = 1, K = .094, and in Dirac space 



(83) 



/(995.,.1) (.0,.0) (.0,.0) (.1,.0) \ 

(.0,.0) (995.,. 1) (.1,.0) (.0,.0) 

[-.2,.0) -(.1,.3) (995., -.1) (.0,.0) 

V(-.l,.3) (.2,.0) (.0,.0) (995., -.1)7 



(84) 



with K = .166, 
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2K 
103 



/ (902., .9) -(1.1, .6) (-.4,.0) (-.9,2.1) \ 

(2.8, .9) (903., 4.2) -(.9,2.1) (-1.9, .0) 

(-2.7, .0) -(1.1,4.1) (902., -.9) (2.8, -.9) 

V(-l.l,4.1) (.7,.0) (-1.1,.6) (903., -4.2)/ 



(85) 



2K 

103 



(11., 5.1) (.04,8.4) \ 



(86) 



In the color off-diagonal sector with K = .166, 

/ -(8.0, 3.6) (5.7,2.1) 
(2.6,-2.3) -(9.8,7.4) -(20., 8.4) -(5.3,8.2) 
-(2.1,7.0) (6.5,1.1) (4.0, .00) (7.3,5.1) 
V (5.9,7.6) (1.4,11.) (3.9,-1.2) (1.3,-13.)/ 
In general, (J)ffuctuations grow as K increases; (J/)in the color singlet sector, scalar com- 
ponents dominate; (J J/) in the color nonsinglet sector, the nonscalar ffuctuations are greater 
than the scalar, sometimes by as much as an order of magnitude; (/V")these statements are 
valid for all A and axial gauge configurations with typical variation between configurations 
of ~ 20%. 

The Euclidean Dirac matrices we use obey Fierz relations 



-pci/9-p7(5 
J- 1 J- 2 



irir2 J- 3 J- 4 

r3r4 



with 



ft>ii fSS ^55" 

Jss — J PP — A ' Jvv 



1 
4 



fSS -I ft!S Q 

Jaa — ^ 5 Jtt — -J 



cSS 



(87a) 



(87b) 



[7^(1 ± ^%Ar{l ± 7')].^ = -[7'^(1 ± 7')]..[7"(1 ± 7')].,- 



If ty = —tx = — t. Wick expansion and Eqs. (|76l), (|T]) and ( gTaHST^) imply 
J = J2{V,s, J(0) dyV^) = Sit) 75 Dt(-t) 7^ , 



x,y 



~ tiix^'L) = riFa = LL- 

■^D72]v/rir2 J'^'^i'^oi'o) otherwise; 



(87c) 



(88) 



(89) 



x,y 







otherwise. 



(90) 
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Comparing (|^) and (|90D to (|88D yields 



(V'oV'o) TiFa = 55; 

"F-,^ = (91) 
^ ^ otherwise; 



r LL; 
CiHl\ = { — (92) 

^^''^ [ -^d72n (V^o^o) /riFa otherwise. 
Therefore, ajj^i^ is the dominant contribution to ao^iiiriVi] unless rir2 = SS, in which 
case affii is 1/N suppressed relative to apiQiss]- By (|87b|) the P = ratios are 



Q;r799 : ar799 : aui'' '■ aw'?'? : ao-w = 1:1:4: —4 : 12 , (93) 

"ss "pp "vv "aa ^tt ' \ ' 

which as described in Ref. are approximately numerically reproduced at /? = 5.7 and 6.0. 
Furthermore, appropriate linear combinations of the give Eq. (|^). Since (ipoipo) is gauge 
invariant, ao±[rir2] is gauge invariant. 
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